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1. Introduction 

Since  December  1971,  the  Aeroelastic  and  Structures 
Research  Laboratory  of  the  Massachusetts  Institute  of  Tech- 
nology has  conducted  research  studies  under  the  sponsorship 
of  the  Air  Force  Office  of  Scientific  Research  on  static  and 
dynamic  nonlinear  structural  problems.  The  general  thread 
of  this  research  effort  is  the  use  of  hybrid  finite  element 
models  for  analyzing  structural  mechanics  problems.  In  addi- 
tion to  the  study  of  problems  with  geometrical  and  material 
nonlinearities,  the  present  task  also  includes  basic  investi- 
gations of  rationalization  and  extensions  of  variational 
formulations  of  finite  element  methods  for  linear  and  nonlinear 
problems.  The  program  also  contains  the  development  of  spe- 
cial crack  elements  for  linear  fracture  mechanics.  Such  ele- 
ments for  both  2-D  and  3-D  problems  are  based  on  the  assumed 
stress  hybrid  model  and  are  derived  by  taking  into  account 
the  singularity  at  the  tip  of  the  crack. 

The  majority  of  results  obtained  under  this  research 
program  have  already  been  documented  either  as  interim  reports 
published  as  AFOSR  Technical  Reports  or  in  the  form  of  techni- 
cal papers  published  either  in  archival  journals  or  in  pro- 
ceedings of  technical  conferences.  The  most  recently  obtained 
results  of  the  following  three  subject  matters  have  not  been 
published: 


(A)  Creep  and  Viscoplastic  Analysis  by  Assumed  Stress 


; 


Hybrid  Finite  Elements 

(B)  Three-Dimensional  Crack  Element  by  Assumed  Stress 
Model 

(C)  Improvement  of  Plate  and  Shell  Finite  Elements  by 
Hybrid  Formulations 

The  first  two  items  are  presented  respectively  as 
Appendices  (A)  and  (B)  of  this  report.  Each  of  these  is  made 
self-contained  with  individual  numbering  systems  for  equations, 

! 

figures  and  references.  The  results  of  the  last  item  has  been 
written  in  a paper  [P-20]* to  be  presented  at  the  AIAA/ASME 
18th  Structures,  Structural  Dynamics  and  Materials  Conference 
at  San  Diego,  California,  March  21-23,  1977. 

In  addition,  the  present  research  program  has  supported 
the  works  of  one  Ph.D.  thesis  and  two  M.S.  theses  in  full 
one  Ph.D.  thesis  in  part  and  one  Ph.D.  thesis  in  progress. 

Also  two  students  completed  their  M.S.  theses  and  one  is 
working  on  his  Ph.D.  thesis  with  some  financial  support  for 
their  computations  under  this  research  program. 

Section  2 is  to  list  all  these  publications,  which  have 
either  already  appeared  or  will  soon  be  in  print. 
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2.  List  of  Publications  under  the  Present  Program 

2.1  Technical  Reports  Published 

1.  Luk,  Chih-Hung,  "Assumed  Stress  Hybrid  Finite  Element 
Method  for  Fracture  Mechanics  and  Elastic-plastic  Analysis", 
AFOSR  TR  73-0493,  M.I.T.  ASRL  TR  170-1,  December  1972,  (also 
M.I.T.,  Department  of  Aeronautics  and  Astronautics,  Ph.D. 
thesis) . 

2.  Lasry , S.J.,  "Derivation  of  Crack  Element-Stiffness  Matrix 
by  the  Complex  Variable  Approach",  AFOSR  TR  73-1602,  M.I.T. 

ASRL  TR  170-2,  February  1973,  (also  M.I.T.,  Department  of 
Aeronautics  and  Astronautics,  M.S.  thesis). 

3.  Lee,  Sung  Won,  "An  Assumed  Stress  Hybrid  Finite  Element 
for  Three-Dimensional  Elastic  Structural  Analysis",  AFOSR  TR 
75-0087,  M.I.T.  ASRL  TR  170-3,  May  1974,  (also  M.I.T.,  Depart- 
ment of  Aeronautics  and  Astronautics,  M.S.  thesis). 

4.  Boland,  P.L.,  "Large  Deflection  Analysis  of  Thin  Elastic 
Structures  by  the  Assumed  Stress  Hybrid  Finite  Element  Method" , 
AFOSR  TR  76-1111,  M.I.T.  ASRL  TR  170-4,  October  1975,  (also 
M.I.T.,  Department  of  Aeronautics  and  Astronautics,  Ph.D. 
thesis) . 

2.2  Technical  Papers  Published 

5.  Pian,  T.H.H.,  "Finite  Element  Methods  by  Variational 
Principles  with  Relaxed  Continuity  Requirement" , Variational 
Method  in  Engineering,  edited  by  C.A.  Brebbia  and  H.  Tottenham 
Southampton  Uni.  Press,  1972,  pp.  3.1-3.24. 
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12.  Pian,  T.H.H.,  "Nonlinear  Creep  Analysis  by  Assumed  Stress 
Finite  Element  Methods",  AIAA  J. , Vol.  12,  No.  12,  December 
1974,  pp.  1756-1758. 

13.  Tong,  P.  and  Pian,  T.H.H. , "Application  of  Finite  Element 
Method  to  Mixed-mode  Fracture",  Recent  Advances  in  Engineering 
Science,  Vol.  6,  1975,  pp.  255-263. 

14.  Pian,  T.H.H. , Spilker,  R.L.  and  Lee,  S.W.,  "Elastic- 
plastic  and  Creep  Analyses  by  Assumed  Stress  Finite  Elements", 
Trans.  3rd,  Int.  Conf.  on  Structural  Mechanics  in  Reactor 
Technology,  Vol.  5,  Part  M,  paper  M 2/1,  pp.  1-7. 

15.  Pian,  T.H.H.  and  Lee,  S.W.,  "Notes  on  Finite  Elements  for 
Nearly  Incompressible  Materials",  AIAA  J. , Vol.  14,  No.  6, 

June  1976,  pp.  824-826. 

16.  Pian,  T.H.H.,  "Hybrid  Models  for  Three-Dimensional  Crack 
Elements",  Proceedings  of  a Workshop  on  Three-Dimensional 
Fracture  Analysis,  edited  by  L.E.  Hulbert,  Battell  Columbus 
Laboratories,  November  1976,  pp.  84,85. 

2.3  Technical  Papers  to  be  Published 

17.  Pian,  T. H.H. , "Variational  Principles  for  Incremental 
Finite  Element  Methods",  to  be  published  in  J.  Franklin 
Institute . 

18.  Boland,  P.L.  and  Pian,  T.H.H.,  "Large  Deflection  Analysis 
of  Thin  Elastic  Structures  by  the  Assumed  Stress  Hybrid  Finite 
Element  Method",  Presented  at  2nd  National  Symposium  on 
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Computerized  Structural  Analysis  and  Design,  Washington,  D.C 
March  29-31,  1976,  (to  be  published  in  Computers  and  Struc- 
tures) . 

19.  Pian,  T.H.H.  and  Boland,  P.L.,  "Formulation  of  Large 
Deflection  Shell  Analysis  by  Assumed  Stress  Finite  Element 
Method",  presented  at  U.S. -German  Symposium  on  Formulations 
and  Computational  Algorithms  in  Finite  Element  Analysis, 
M.I.T.,  August  9-13,  1976,  (to  be  published  by  M.I.T.  Press) 

20.  Lee.  S.W.,  and  Pian,  T.H.H. , "Improvement  of  Plate  and 
Shell  Finite  Elements  by  Hybrid  Formulations",  AIAA  paper 
No.  77-413.  For  presentation  at  AIAA/ASME  18th  Structures, 
Structural  Dynamics  and  Materials  Conference,  San  Diego,  CA, 
March  21-23,  1977. 

2.4  Theses  Supported  in  Part  in  the  Present  Program 

21.  Shimizu,  T. , "Axisymmetric  Creep  Analysis  by  Assumed 
Stress  Hybrid  Finite  Element  Method",  M.I.T. , Department  of 
Aeronautics  and  Astronautics,  M.S.  thesis,  September  1974. 

22.  Scharnhorst,  Thomas,  "Variational  and  Finite  Element 
Formulations  of  Rubber-Like  Materials",  M.I.T.,  Department 
of  Aeronautics  and  Astronautics,  M.S.  thesis,  June  1976. 
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3.  Summary  of  Research  Findings 


The  research  results  obtained  under  this  program  are 
classified  into  the  following  four  categories: 

(A)  Rationalization  and  extensions  of  variational  formu- 
lations of  finite  element  methods 

(B)  Development  of  crack  elements 

(C)  Analysis  of  geometrically  nonlinear  problems 

(D)  Analysis  of  material  nonlinear  problems 

The  research  findings  under  these  categories  are  summarized 
in  the  following. 

3.1  Rationalization  and  Extensions  of  Variational  Formu- 
lations of  Finite  Element  Methods 

(A)  In  the  application  to  finite  element  formulations  in 
solid  mechanics  it  has  been  demonstrated  [P-5]  that  the  con- 
ventional variational  principles  can  be  modified  and  extended 
by  allowing  the  relaxation  of  the  continuity  requirements 
along  the  interelement  boundary  and  that,  correspondingly, 
many  versions  of  hybrid  models  can  be  constructed.  One 
example  is  that  the  so-called  assumed  stress  hybrid  elements 
can  be  derived  using  either  a modified  complementary  energy 
principle  [R-l,  R-2]'*'  or  the  Hellinger-Reissner  principle  if 
in  the  latter  the  assumed  stresses  are  also  in  equilibrium  [P-5]. 


^Reference  numbers  headed  by  R refer  to  that  in  the  reference 
list  at  the  end  of  this  report. 
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The  latter  approach  is  computationally  more  efficient,  parti- 
cularly for  the  derivation  of  3-d  elements  [P-3].  An  eight- 
node  isoparametric  3-D  assumed  stress  element  has  been  found 
to  need  a comparable  amount  of  computing  time  for  its  genera- 
tion in  comparison  with  the  conventional  assumed  displacement 
model,  but  to  provide  more  accurate  results.  Another  example 
is  the  utilization  of  a modified  Hellinger-Reissner ' s principle 
to  derive  the  geometric  stiffness  matrix  for  buckling  analysis 
and  the  mass  matrix  for  vibration  analysis  for  assumed  stress 
hybrid  elements  [P— 10 ] . 

(B)  Variational  principles  and  modified  principles  in 
solid  mechanics  are  extended  to  cover  incremental  finite 
element  methods  [P-17].  For  those  principles  based  on  com- 
plementary energy,  the  assumed  stresses  need  not  satisfy  the 
complete  equilibrium  conditions.  Corrections  for  effects  of 
equilibrium  imbalance  and  compatibility  mismatch  at  the  be- 
ginning of  each  loading  increment  can  be  taken  into  account 
in  a systematic  manner. 

(C)  It  has  been  found  that  hybrid  formulations  can  be 
developed  to  reduce  severe  constraints  that  appear  in  the 
derivation  of  several  finite  elements  by  the  conventional 
assumed  displacement  method,  hence  can  lead  to  more  satisfac- 
tory elements.  In  the  assumed  stress  hybrid  model,  for 
example,  the  constraint  for  incompressibility  is  reduced  to 
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only  one  condition,  hence  an  ordinary  hybrid  stress  element 
can  be  applied  to  problems  which  are  very  near  to  completely 
incompressible  [P-15].  Two  other  examples  are  the  use  of 
hybrid  elements  to  extend  the  range  of  applications  for  plate 
and  shell  problems  [P-20]. 

3.2  Development  of  Crack  Elements 

A study  has  shown  that  the  convergence  rate  of  the  finite 
element  method  for  problems  with  singularities  is  very  slow 
if  ordinary  elements  are  used  [P-6].  Thus,  it  is  desirable 
to  employ  special  crack  elements  in  which  the  proper  singu- 
larities are  taken  into  consideration.  Two  basic  types  of 
assumed  stress  hybrid  models  have  been  used  to  derive  crack 
elements  for  plane  elasticity  problems.  In  the  first  one, 
the  assumed  stresses  contain  the  correct  distribution  of  the 
singular  term  but  the  additional  terms  consist  of  only 
polynomial  expansions  in  the  radial  distance  r from  the  crack 
tip  [ P— 1 ] . In  this  case,  several  crack  elements  are  needed 
around  the  crack  tip.  In  the  second  case,  the  assumed  stresses 
consist  of  complete  series  expansion  of  the  exact  solution 
around  the  crack  tip.  Thus,  all  terms  satisfy  both  stress 
equilibrium  and  compatibility  conditions  and  that  contain 
not  only  the  singular  terms  and  polynomial  expansion,  but  also 
terms  involving  rn ^ with  n being  integers.  In  this  case, 
the  special  element  is  one  single  element  with  an  inbedded 
crack.  [P-2,  P-7,  P-8,  P-13]. 
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Of  these  two  types  of  crack  elements,  the  second  type  is 
apparently  more  efficient  both  in  the  construction  of  the 
element  and  in  the  accuracy  of  the  solution.  Both  elements 
have  been  demonstrated  to  be  superior  among  the  existing 
finite  element  methods  for  the  determination  of  stress  in- 
tensity factors  [ R—  3 ] . The  second  approach  for  hybrid  crack 
elements  has  been  extended  under  other  ASRL  research  programs 
sponsored  by  the  Air  Force,  to  problems  involving  bi-material 
interface  [R-4,  R-5].  A nine-node  hybrid  crack  element 
with  a code  name  of  PCRK  59  has  been  employed  to  compute 
stress  intensity  factors  for  several  aircraft  structural 
details  including  attachment  lugs,  centered  and  offset  fastener 
holes  [R-6 , R-7,  R-8,  R-9].  This  element  is  now  incorporated 
in  the  library  of  subroutines  of  the  FEABL-".  software  [R-10] 
and  has  also  been  incorporated  by  several  aerospace  manu- 
facturers in  their  structural  analysis  computing  programs. 

The  first  approach  for  hybrid  crack  elements  has  been 
extended  to  3-D  crack  analysis  [P-16].  A brief  summary  of 
this  development  is  given  in  Appendix  B of  this  report. 

3.3  Analysis  of  Geometrically  Nonlinear  Problems 
(A)  The  bifurcation  and  the  postbuckling  behaviors  of 
shell  of  revolution  have  been  examined  by  the  finite  element 
method  [P-9].  The  symmetric  prebuckling  deformation  of  the 
shell  is  first  formulated  and  the  stiffness  matrices  of 
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various  harmonics  are  examined  for  the  possibility  of  bifurca- 
tion. A postbuckling  analysis  is  then  performed  to  examine 
the  stability  of  such  bifurcated  mode.  The  effect  of  symmetric 
and  asymmetric  imperfections  and  of  material  orthotropy  on 
the  buckling  load  of  spherical  shells  has  been  studied. 

(B)  The  assumed  stress  hybrid  finite  element  model  has 
been  demonstrated  to  be  suitable  for  incremental  solution  of 
thin  elastic  structures  in  large  deflections.  An  incon- 
sistent model  may  be  adopted  in  place  of  a consistent  one 
which  would  involve  a much  more  complicated  task  of  assuming 
stresses  to  satisfy  the  entire  equilibrium  equations.  The 
discrepancy,  however,  can  be  taken  into  account  by  the  inclu- 
sion of  correction  terms  corresponding  to  checks  on  stress 
equilibrium  and  compatibility  at  the  beginning  of  each  loading 
increment.  The  method  has  been  applied  to  plate  and  shells 
in  large  deflections  but  in  small  strains  [P-4,  P-18,  P-19] 
and  to  rubber-like  material  in  the  finite  strain  range  (P-22]. 

3.4  Analysis  of  Material  Nonlinear  Problems 

(A)  For  elastic-plastic  analyses,  it  has  been  shown  that 
both  the  tangent  stiffness  approach  and  the  initial  stress 
approach  can  be  adapted  in  conjunction  with  the  assumed  stress 
hybrid  model  [P-1,  P-11,  P-14].  Results  have  demonstrated  that 
the  hybrid  elements  are  more  efficient  than  the  conventional 
assumed  displacement  element  from  the  point  of  view  of  solu- 
tion accuracy  and  computing  effort. 
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(B)  The  assumed  stress  hybrid  elements  have  also  been 
applied  to  nonlinear  creep  analysis  (P-12,  P-14,  P-21].  The 
resulting  equations  are  in  the  form  of  a system  of  nonlinear 
first  order  deferential  equations  for  which  many  standard 
numerical  techniques  are  available.  Appendix  A presents 
some  results  and  also  demonstrates  that  the  elastic-plastic 
problems  can  be  solved  by  the  viscoplasticity  approach,  hence 
can  also  be  solved  by  the  present  creep  analysis  program. 
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APPENDIX  A 

i 

Creep  and  Viscoplastic  Analyses  by  Assumed 
t Stress  Hybrid  Finite  Elements 

T.  H.  H.  Pian  and  S.  W.  Lee 
ABSTRACT 

A hybrid  stress  finite  element  formulation  is  derived 
for  creep  analysis  and  elastoplastic  analysis  by  viscoplastic 
formulation.  Application  of  the  Hellinger-Reissner  principle 
leads  to  a system  of  nonlinear  first  order  ordinary  differen- 
tial equations  with  stress  parameters  in  each  element  as 
variables.  Examples  involving  four-node  plane  stress  element 
and  solid  of  revolution  element  indicate  that,  for  elements 
with  a low  order  displacement  field,  the  assumed  stress  hybrid 
formulation  is  computationally  more  efficient  than  the  con- 
ventional assumed  displacement  formulation. 

I.  Introduction 

* Creep  and  viscoplastic  problems  are  nonlinear  in  nature 

and  in  general  require  numerical  solution  methods.  The  finite 
element  method  based  on  the  assumed  displacement  models  has 
been  applied  to  this  type  of  problem  treating  creep  or  visco- 

} 

plastic  strains  as  initial  strains.  Percy  et  al  [1], 

I 


Greenbaum  and  Rubinstein  [2],  Sutherland  [3],  Chang  and  Rashid 
[4],  Branca  and  Boresi  [5],  Donea  and  Giuliani  [6],  etc. 
solved  small  displacement  creep  problems  while  Cyr  and  Teter 
[7]  solved  large  deformation  creep  problems  combined  with 
plastic  effect.  Numerical  stability  is  very  important  for  all 
numerical  analyses  and  sufficiently  small  time  steps  were 
used  in  these  works.  Cormeau  [8]  developed  a formulation  for 
viscoplastic  analysis  of  solids  utilizing  the  constitutive 
law  proposed  by  Perzyna  [9]  and  applied  it  to  elastoplastic 
analysis  and  to  analysis  of  transient  creep  under  the  time 
hardening  law.  A theoretical  stability  criterion  has  been 
established  in  his  formulation.  Bodner  et  al  [10,11]  pro- 
posed an  alternative  viscoplastic  constitutive  law  and  applied 
it  to  the  finite  element  method. 

In  this  work  an  assumed  stress  hybrid  finite  element 
method  for  small  deformation  creep  analyses  will  be  derived 
from  the  Hellinger-Reissner  principle  with  initial  strains. 

The  resulting  system  of  nonlinear  first  order  ordinary  dif- 
ferential equations  can  be  solved  by  various  Runge-Kutta 
methods.  Finally  by  employing  Perzyna ' s constitutive  law  for 
viscoplasticity,  the  timewise  solution  procedure  is  to  be 
used  for  solution  of  elastoplastic  problems. 


16 


✓ 


!.  Constitutive  Equations 
(a)  Creep  Law 

For  the  present  purpose,  phenomenological  description 
of  creep  will  be  adopted.  In  the  study  of  long  duration 
creep  the  steady  state  creep  law  is  used  and  the  corresponding 
uniaxial  creep  strain  rate  is  expressed  by 

i C = F C S',  T)  (1) 

where  a is  the  uniaxial  stress  and  T is  the  temperature. 

If  the  effects  of  a and  T are  assumed  to  be  separable, 

then 

ec  = -f CO  ^Ct)  (2) 


For  the  stress-dependent  function,  we  use  Norton's  power 
law , i . e . , 


-KO  = &6~'"  <3> 

where  B and  n are  material  constants.  For  the  present  study 
it  will  be  assumed  from  now  on  that  the  temperature  remains 
constant.  In  the  case  of  transient  creep,  the  time  hardening 
law  and  the  strain  hardening  law  are  available.  Under  constant 


17 


r 


I 


r 


t 


1 


temperature  and  constant  stress  the  uniaxial  creep  strain 
may  be  expressed  as 

£ = A <5“"  t*  (4) 

where  A,  m,  k are  material  constants  and  t represents  time. 
Differentiating  eq.  (4)  with  respect  to  time  leads  to 

£c  = A-*.  (5) 

The  variable  t can  be  eliminated  from  eq.  (5)  by  solving  eq. 
(4)  for  t and  substituting  into  eq.  (5).  The  result  is 

ic='k6'J*  A*  (£e)’  V 

Equations  (5)  and  (6)  are  called  the  time  hardening  law  and 

the  strain  hardening  law  respectively.  The  time  hardening 

law  can  be  written  in  a form  which  resembles  the  steady  state 

k , 

creep  law  by  defining  a parameter  x = t and  differentiating 
eq.  (4)  with  respect  to  x to  yield 

r = as"*  (7) 


Under  multiaxial  stress  condition  creep  strain  rates 
are  expressed  as 


f 


18 


where  eC  is  the  equivalent  creep  strain  rate  and  the  yield 
e 


condition  is  represented  by 


FC  ° (9) 

By  substituting  the  uniaxial  creep  strain  rates  of  eqs . (6) 

and  (7)  into  eq.  (8)  for  the  equivalent  creep  strain  rate  and 
using  the  equivalent  stress  oe  in  place  of  the  uniaxial 
stress,  we  obtain 


dF 


(10) 


for  the  steady  state  creep  or  creep  under  the  time  hardening 
law  and 


A 

= tec+ 


A*ce;>'  'S : 


i - 


i 


(ID 


for  creep  under  the  strain  hardening  law.  In  eq.  (10)  the 
creep  strain  rate  is  determined  by  the  present  stress  values 
while  in  eq.  (11)  it  is  expressed  in  terms  of  the  present 
stresses  and  the  equivalent  creep  strain  or  creep  strains. 


If  the  Mises-Hencky  yield  criterion  is  used,  then 


3F  , 3 5tj 


(12) 


where  are  the  deviatoric  stresses. 

(b)  Perzyna's  Viscoplastic  Law 

For  materials  undergoing  viscoplastic  process,  multi- 


axial  viscoplastic  strain  rates  are  expressed  as  follows 


H. 


(13) 


where  F is  the  same  yield  function  as  in  eq.  (9)  with  uni- 
axial yield  stress  and 


<$>  = 6*c  “ Sy 

if 

IV 

5V 

= o 

if 

< 

(14) 


and  y is  a material  constant.  If  the  yield  stress  is  constant 
the  viscoplastic  strain  rates  are  determined  by  stresses  only. 
The  above  representation  has  been  proposed  by  Perzyna  t9). 

It  is  different  from  the  classical  inviscid  plasticity  theory 
in  that  it  permits  the  stresses  to  exceed  the  yield  stress. 

The  visco-plastic  formulation  can  be  used  to  solve  classi- 
cal plasticity  problems.  This  is  achieved  by 


* 


finding  a steady  state  solution  under  a given  load  incre- 
ment. 

(c)  Bodner-Partom' s Law 

In  Bodner-Partom ' s constitutive  law  inelastic  strain  rates 
are  expressed  as  follows 


| 


£v/  (15) 

where  is  the  second  invariant  of  the  deviatoric  stresses 
and  the  function  f include  material  constants  and  plastic 
work  done.  No  yield  stress  is  defined  and  thus  inelastic 
strain  components  are  present  all  the  time  from  the  moment 
a load  is  applied.  In  the  present  investigation  we  will  not 
use  this  constitutive  law.  But  it  can  be  treated  in  the 
same  manner  as  creep  laws  or  Perzyna's  law. 


3.  Finite  Element  Formulation 

An  assumed  stress  hybrid  finite  element  can  be  formulated 
by  a modified  complementary  energy  principle  or  by  the 
Hellinger-Reissner  principle.  For  the  example  problems 
carried  out  in  this  investigation,  the  Hellinger-Reissner 
principle  is  suitable  since  it  is  easy  to  construct  compatible 
displacement  fields.  For  creep  and  viscoplastic  analyses, 
the  Hellinger-Reissner  principle  is  extended  to  include  the 


21 


i 


initial  strains.  The  functional  it  is  written  in  terms  of 

JR 

stress  and  displacement  rates  as  follows. 

- 5 j [ <tvj  i ( itj  + ij.i ) - i 

- -PiUtJJv  - / (] 


where 


o.  . = stress  rate  tensor 

uk  = displacement  rates 

S . . , „ = elastic  coefficient  tensor 
l^kJl 

= initial  strain  rates,  i.e.,  creep  strain  or 
viscoplastic  strain  rates  in  the  present  case 

F.  = applied  body  force  rates  per  unit  volume 

• 

T^  = applied  traction  rates 

V = volume  of  the  element 
n 

S = surface  of  the  nth  element  over  which  tractions 

o 

n 

are  applied 

In  the  matrix  form  this  functional  can  be  written  as 

- i iTl  £ - £i' 

” C 

- Ftk] ^ - $ t"ids 

The  above  expression  is  valid  for  small  deformation.  The 
corresponding  variational  principle  for  large  deformation 
has  been  given  by  Sanders  et  al  [12]  and  Pian  [13].  In 
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eq.  (16)  the  displacements  and  the  stresses  are  independent 
quantities.  For  the  hybrid  stress  finite  element  implementa- 
tion, stress  and  dispalcement  rates  are  assumed  in  terms  of 
unknown  stress  rate  parameters  8 and  nodal  dispalcement  rates 
q respectively,  i.e.. 


S'  = ? /S  + Pf 

A/  ^ <v' 

IX  =.  A hr 

/v  ^ 

(18) 

then  (£!l)=£i 


Here  the  assumed  stress  rates  are  made  to  satisfy  the 
equilibrium  condition.  Thus  P 8 is  the  homogeneous  solution 
of  the  equilibrium  equations  and  P„  8_  is  a particular  solu- 

~ r ~ r 

tion  of  the  equilibrium  equations.  Substituting  eq.  (18) 
into  eq.  (17)  and  performing  necessary  integrations,  we  obtain 

TTR  «z  ( i + <LT  (tf  % - i li 

~ /iT  Hf  &-  - |3T  ) 

= ( pFT  0 Jv 

A*  -'V  A*  A* 

•W 

= $„  tTS  Pf'1'' 

iTe.  = 5 ttu  <is 


<T  = J PT  B d V 

^ ^ (I  /v» 

y> 

< = C PT  ll’dv 


(20) 


Taking  6 ttr 


0 with  respect  to  3 gives 
<k%-  - H fl  - Hpi  - i = o 


or 

8 = H'1  C St  - - S-l 

rv  v 

Substituting  eq.  (21)  into  eq.  (19)  yields 


(21) 


iT 


cs. 


(22) 


where 


GTH  1G  = element  stiffness  matrix 

9a  + (9T?~1?f  ' 5fT)?f 


(23) 


= rate  of  applied  nodal  load 

6°  = GTH-'*'G„  = rate  of  equivalent  nodal  load  due  to 
~n  ~ ~ ~0 

inelastic  strain 

The  element  stiffness  matrix  and  nodal  load  can  be  assembled 
to  yield 


. T r * 

£ l £ 


■+ 


(24) 


and  the  stationariry  condition  of  irR  with  respect  to 


q gives 


k i = & + 


(25) 
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Solving  eq.  (25)  for  q and  substituting  into  eq.  (21) , we 
obtain 


i = H-  ( S is'gT  H-  -I)  j.  - tCt  !£*&  - £'  Sp  b 

For  problems  involving  state  creep,  creep  under  the 
dening  law  and  viscoplasitc  law  with  constant  yield 

V '±^  = ±(6) 


and  thus  from  eq.  (26) 

i - i(sl 


(26) 

time  har- 
stress . 

(27) 

(28) 


On  the  other  hand,  for  creep  under  the  strain  hardening  law 
or  for  viscoplasticity  law  with  yield  stress  dependent  on 
the  viscoplastic  strains 


r 


Equations  (28)  and  (29)  are  systems  of  nonlinear  first  order 
ordinary  differential  equations  with  8 and  8,e°  as  unknowns 
respectively.  In  the  finite  element  formulation,  integrations 
over  an  element  are  performed  numerically.  Thus,  in  eq.  (29) 
the  vector  e°  represents  inelastic  strains  at  Gaussian  inte- 
grations points.  The  assumed  stress  hybrid  formulation  for 
the  steady  state  creep  and  for  transient  creep  under  the  time 
hardening  law  has  been  given  by  Pian  [14]. 

In  the  conventional  assumed  displacement  method  based  on 
the  principle  of  virtual  work  with  initial  strains,  the  un- 
knowns are  stresses  a at  Gaussian  integration  points  instead 
of  8 and  we  have  the  same  type  of  equations  as  eqs.  (28)  and 
(29)  with  8 replaced  by  a.  The  amount  of  computing  time 
depends  on  the  number  of  unknowns.  For  elements  with  a low 
order  assumed  displacement  field,  assumed  stress  hybrid  for- 
mulation will  result  in  a problem  with  fewer  number  of  un- 
knowns than  that  by  the  corresponding  assumed  displacement 
method.  For  example,  consider  a four  node  plane  stress 
element  derived  by  2 x 2 integration  rule.  An  element  based 
on  the  displacement  approach  will  have  three  stress  components 
at  each  Gaussian  point  and  thus  has  twelve  unknown  stresses 
per  element.  On  the  other  hand  a hybrid  stress  element 
with  7 8's  has  seven  unknowns.  Similarly,  for  the  assumed 
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stress  hybrid  model,  a four-node  axisyrrunetric  element  with 
8 g's  and  a eight-node  three  dimensional  element  with  24  8's 
have  fewer  number  of  unknowns  that  the  corresponding  ele- 
ments by  the  assumed  displacement  model.  This  advantage  of 
assumed  stress  hybrid  elements  disappear  as  the  number  of  nodes 
are  increased  for  an  element  since  more  g's  are  needed  to 
prevent  hazardous  kinematic  modes.  For  example  for  an  eight 
node  plane  stress  element  the  number  of  unknown  fj's  in  the 
hybrid  formulation  and  the  number  of  o's  for  the  displacement 
formulation  are  comparable  to  each  other. 

4.  Solution  Method 

Equations  (28)  and  (29)  can  be  solved  timewise  by  numeri- 
cal methods  such  as  Runge-Kutta  methods  [15,16]  using  elastic 
solution  as  initial  conditions.  All  these  methods  have  numeri- 
cal stability  limit.  For  the  integration  of  eq.  (28)  stabi- 
lity criterion  developed  by  Cormeau  [8]  has  been  utilized. 

The  well  known  fourth  order  explicit  Runge-Kutta  method  was 
used  for  steady  state  creep  and  creep  under  the  time  hardening 

law.  In  the  case  of  the  time  hardening  law,  the  integration 

k.  • 

is  performed  for  the  time  parameter  x = t . The  real  time  is 
recovered  by  t = x^^.  On  the  other  hand  the  midpoint  Runge- 
Kutta  method  is  considered  more  appropriate  for  elastoplastic 

*~r 

problems  where  a lower  order  method  was  chosen  due  to  the 
possibility  of  piecewise  continuity  of  stresses. 
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To  solve  elastoplastic  problems,  it  is  assumed  at 

first  that  the  present  state  is  on  the  yield  surface  or  in 

an  elastic  state.  Given  a load  increment,  the  material  will 

flow  according  to  eq.  (13)  if  o > o . But  if  the  structure 

6 y 

is  stable,  the  plastic  flow  will  stop  after  a certain  time 
and  a steady  state  is  reached  under  the  given  load.  This 
steady  state  will  be  the  one  which  is  obtainable  by  the 
classical  plasticity  formulation.  In  actual  calculations  it 
is  assumed  that  a steady  state  has  been  reached  if 


6e  ~ 


< £ 


( 30) 


where  e is  an  arbitrarily  assigned  value.  For  the  present 
numerical  solution  it  is  chosen  as  0.01. 

For  creep  problems  under  the  strain  hardening  law,  the 
time  increment  for  the  numerical  integration  of  eq.  (29)  has 
been  determined  such  that 

ec  At  < total  strain  at  time  t x a (31) 

e 


where  a is  a preassigned  value.  The  value  of  a = 0.05  was 
used  in  the  present  work.  It  is  to  be  noted  that  for  materials 
with  k < 1,  eq.  (11)  and  hence  eq.  (29)  are  singular  at 
time  t = 0 since  = 0.  Therefore  it  is  necessary  to  use 
a method  which  does  not  require  evaluation  of  the  creep  strain 
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rate  at  t = 0.  A fourth  order  implicit  Runge-Kutta  method 
developed  by  Butcher  [17]  was  chosen  for  this  purpose.  After 
one  time  increment,  it  is  switched  to  the  fourth  order 
explicit  method  since  the  explicit  method  requires  less  com- 
puting time  for  the  same  accuracy. 


5.  Numerical  Examples 

Two  assumed  stress  hybrid  elements  are  used  for  numerical 
solution  of  example  problems.  They  are  a four-node  quadri- 
lateral plane  stress  element  with  the  following  stress 
assumption 

X + /3S  £ 

3-  + * (32) 

= 6,  - a4  3.  - 6,  * 

and  a four-node  quadrilateral  axisymmetric  element  [191  with 
the  following  assumed  stresses 

6V  * 6,  + + 6i  */r 

<Te  = (33) 

€TZ  = E>  g + 6t  /r  - V-r 

sri  = - B,  + &4  «/r  + «7  + /r 

The  isoparametric  representation  is  used  for  displacement 
assumption.  Among  the  following  examples,  the  first  two 
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problems  deal  with  creep  and  the  last  one  with  elastic- 
plastic  analysis.  Computations  were  performed  with  an  IBM 
370/165  machine  at  the  M.I.T.  Information  Processing  Center. 

(a)  Creep  of  a thick  cylinder  under  internal  pressure. 

An  infinitely  long  thick  cylinder  subjected  to  internal 
pressure  was  solved  as  a plane  strain  problem  using  six  axi- 
symmetric  elements.  The  geometry,  load  and  material  pro- 
perties are  given  in  fig.  1.  Results  are  given  in  figs.  1 
and  2.  The  time  hardening  law  and  the  strain  hardening  law 
gave  almost  identical  results.  All  results  are  close  to 
those  obtained  by  Greenbaum  and  Rubinstein  [2]. 

(b)  Creep  of  a rotating  disk 

This  is  a problem  involving  distributed  body  forces. 

The  geometry  and  material  constants  are  given  in  fig.  3. 

First  it  was  analyzed  by  a row  of  six  axisymmetric  elements  and 
then  by  a row  of  twelve  plane  stress  elements.  In  the  plane 
stress  element  a particular  solution  of  the  equilibrium 
equations  was  chosen  for  the  averages  of  the  body  forces 
over  the  volume  of  the  element.  For  the  axisymmetric  ele- 
ment the  particular  solutions  were  simply  set  to  zero. 
Discussions  relating  these  two  approaches  are  given  in  ref. 

18.  Calculated  results  for  the  time  hardening  law  are  shown 
in  figs.  3 and  4.  In  fig.  3,  t and  m represent  the  thick- 
ness and  the  density  of  the  disk  respectively. 
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(c)  An  infinitely  long  thick  cylinder  of  elastic-perfectly 
plastic  material  under  internal  pressure  p 
The  geometry  and  material  properties  are  as  follows. 


Inner  radius  a 


5 inch. 


Outer  radius  b =10  inch. 

Young's  modulus  E = 10^  psi 

Poisson's  ration  v =0.3 

Yield  stress  a = 20000  psi 

Table  1 shows  the  radial  displacement  at  the  inner  wall  cal- 
culated by  both  the  assumed  stress  hybrid  method  and  the 
assumed  displacement  method  using  eight  axisymmetric  elements. 
The  2x2  Gaussian  integration  points  were  used  for  numerical 
integration.  The  exact  solutions  were  obtained  according  to 
the  procedure  described  in  Ref.  [21]  based  on  the  classical 
plasticity  theory.  The  accuracy  of  the  finite  element  solu- 
tion decreases  as  the  load  approaches  collapse  value  and 
more  fictitious  time  steps  are  needed  to  reach  a steady 
state  under  the  given  load  increment.  But  it  is  seen  that 
the  assumed  stress  hybrid  method  gives  better  solution  in  less 
computing  time  than  the  assumed  displacement  method. 


6.  Conclusion 

An  assumed  stress  hybrid  finite  element  formulation  has 
been  derived  for  viscoplastic  and  creep  analyses.  The  prac- 
ticality of  the  method  has  been  demonstrated  by  example 
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problems.  For  elements  with  a low  order  of  assumed  dis- 
placement field  the  assumed  stress  hybrid  formulation  is 
computationally  more  efficient  than  the  conventional  dis- 
placement method  for  small  deformation  problems. 
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Table  1.  Radial  Displacement  x 10  at  the  Inner  Wall 


p/a 

v y 

Exact 

Solution 

Hybrid 

Method 

Displacement 

Method 

0.5117 

0.1002 

0.1006 

0.0995 

0.5786 

0.1206 

0.1210 

0.1191 

0.6345 

0.1435 

0.1432 

0.1416 

0.6806 

0.1689 

0.1688 

0.1665 

0.7179 

0.1966 

0.1965 

0.1923 

0.7472 

0.2264 

0.2255 

0.2191 

0.7692 

0.2583 

0.2550 

0.2474 

0.7846 

0.2922 

0.2871 

0.2746 

0.7937 

0.3278 

0.3176 

0.2963 

0.7972 

(collapse) 

0.3650 

0.3328 

0.3061 

CPU  Time 

(sec. ) 

6.30 

9.29 

IISPLACEMENT  AT  THE  INNER  WALL  VS.  TIME 


0.16 


0.18 


0.20 


0.22 


0.24  in. 


FIG.  2 STRESS  VS.  RADIUS 


FIG.  3 RADIAL  DISPLACEMENT  AT 


APPENDIX  B 


IT 

v,'\ 


* 

I 


, 


I 


I 

j' 

. 

l l 

4 


i 


I 


Three-Dimensional  Crack  Element  by 
Assumed  Stress  Hybrid  Model 

T.  H.  H.  Pian  and  K.  Moriya 

1.  Introduction 

There  exist  three  basic  finite  element  hybrid  models 
for  the  development  of  special  singular  elements  for  linear 
fracture  mechanics  [1].  The  basic  schemes  for  these  models 
are  as  follows. 

(1)  A scheme  which  is  based  on  the  assumed  equili- 
brating stress  field  which  contains  the  singular  stress  term 
and  also  satisfies  the  compatibility  condition  inside  the 
element,  and  on  independently  assumed  boundary  displacements. 
Such  model  may  be  interpreted  as  either  a hybrid  stress 
model  or  a hybrid  displacement  model. 

(2)  A scheme  which  is  based  on  assumed  equilibrating 
stress  field  including  the  singular  term  and  on  independently 
assumed  boundary  displacements.  Such  model  is  a hybrid 
stress  model. 

(3)  A scheme  which  is  based  on  assumed  displacement 
field  including  the  singular  term  and  on  independently  assumed 


41 


r 


boundary  displacements  and  boundary  tractions.  Such  a 
model  is  a hybrid  displacement  model. 

The  applications  of  these  three  schemes  to  two  dimen- 
sional crack  elements  have  been  reported  by  Tong,  Pian  and 
Lasry  [2],  by  Pian,  Tong  and  Luk  [3],  and  by  Atluri,  Kobayashi 
and  Nakagaki  [4,5]  respectively.  It  appears  that  for  two- 
dimensional  problems  the  most  desirable  element  is  one 
formulated  by  the  first  scheme  for  which  only  one  crack  ele- 
ment is  needed  and  an  imbedded  crack  is  included  in  the  ele- 
ment. In  this  case  only  the  nodal  displacements  are  left  as 
unknowns.  In  implementing  such  element,  only  boundary 
integrations  are  called  for.  Since  the  element  boundary 
does  not  contain  the  crack  tip,  the  integration  does  not  in- 
volve any  singular  terms.  The  application  of  the  other 
two  schemes,  however,  will  result  in  the  use  of  several  ele- 
ments around  the  tip  of  the  crack.  It  is  then  needed  to  intro- 
duce the  stress  intensity  factors  as  unknowns  in  addition  to 
the  nodal  displacements.  Furthermore,  the  evaluation  of 
the  stiffness  matrix  in  these  cases  involves  boundary  inte- 
gration with  singular  terms  and  special  treatments  are 
needed  in  the  numerical  integration  process.  Of  the  last  two 
approaches,  the  hybrid  stress  model  involves  only  two  sets 
of  field  variables,  while  the  hybrid  displacement  model  con- 
tains three  sets  of  field  variables. 
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For  three-dimensional  fracture  mechanics,  it  is  unfor- 
tunate that  the  first  approach  is  no  longer  applicable.  First 
of  all  the  only  available  near  field  solution  near  a crack 


front  is  a two-dimensional  behavior  in  terms  of  r and  0, 
where,  as  shown  in  fig.  1,  r is  the  radial  distance  from  the 
crack  front,  and  0 is  the  angle  from  the  n-t  plane.  Here  the 
t-axis  is  tangent  to  the  crack  front,  the  n-axis  is  on  the 


plane  of  the  crack  and  z is  normal  to  the  crack  plane. 

The  singular  behavior  of  the  stress  field  and  the  corres- 
ponding displacement  distribution  are  given  by  the  following 
equations . 
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Here,  K^,  and  are  the  stress  intensity  factors  and 
are  of  unknown  distributions  along  the  crack  front.  Thus, 
a complete  series  solution  with  each  term  as  known  function  of 
r,  0 and  t is  not  readily  available.  Also,  since  the  crack 
front  will  always  intersect  the  boundary  of  the  element,  it  is 
necessary  in  the  formulation  of  the  element  stiffness  matrix 
to  evaluate  some  surface  integrals  that  are  singular.  Atluri 
et  al  have  applied  the  third  approach,  or  the  hybrid  displace- 
ment model  to  derive  three-dimensional  crack  elements.  The 
present  note  is  to  describe  the  development  of  three-dimensional 
crack  elements  using  the  second  approach  or  the  hybrid  stress 
model. 

2.  Three-Dimensional  Crack  Elements  by  Assumed  Stress  Hybrid 
i Model 

For  the  presently  derived  three-dimensional  elements 
the  edges  are  all  made  straight  lines.  Thus,  if  in  a problem 
the  crack  front  is  curved  it  must  be  approximated  by  an  assem- 
blage  of  straight  segments.  Each  segment  of  the  crack 
front,  then,  serves  as  the  common  edge  of  a group  of  brick 
type  special  crack  front  elements.  An  example  arrangement  with 
four  such  elements  is  shown  in  Fig.  2.  These  elements  are 
classified  as  type  A element  which  contains  the  crack  surface 
and  type  B element  which  does  not  contain  the  crack  surface. 


elements  labelled  as  A'  and  B'  are  those  with  the  crack  front 
located  at  the  upper  corners  of  the  block  elements.  In 
analyzing  a complete  solid  with  general  mesh  pattern,  regular 
finite  elements  are  used  to  surround  these  special  elements. 

The  variational  functional  ir  to  be  used  in  the  derivation 
of  this  element  is  written  as 


n(aij'  ui}  n=|;Vn2  Cijkl  aij  °kl 
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where  the  stress  ck  ^ satisfies  equations  of  equilibrium, 

C.  „ is  the  elastic  compliance  tensor,  T.  is  the  surface 
ijk£  1 

traction  obtained  from  the  stresses  o^,  u^  is  the  boundary 
displacement  and  is  prescribed  surface  traction  over  S„  . 
Thus,  the  independent  field  variables  are  the  equilibrating 
stresses  and  boundary  displacements  CL.  In  a crack  front 

element,  the  assumed  stress  includes  not  only  the  regular 
polynomial  terms  but  also  the  asymptotically  exact  term  which 
is  derived  from  the  well  known  local  solution  of  an  embedded 
elliptical  crack,  and  is  given  in  eq.  (1).  Hence,  the 
assumed  stress  in  the  crack  fron  element  is  written  as 


The  first  term  of  the  above  equation  corresponds  to 
the  singular  solution  eq.  (1) , while  the  second  term  contains 
only  regular  polynomials.  Both  terms  satisfy  homogeneous 
equilibrium  equations  and  the  stress  free  conditions  over  the 
crack  surfaces. 

The  assumed  boundary  displacement  CL  which  is  expressed 
in  terms  of  nodal  displacement  contains  asymptotically  correct 
variation  of  displacement  on  the  faces  that  meet  or  intersect 
the  crack  front  line,  and  the  condition  of  displacement  com- 
patibility is  completely  satisfied  across  all  interelement 
boundaries.  Thus,  CL  is  expressed  as 


a. 
u . 
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where  q^  is  a generalized  (nodal)  displacement  and  ^ is  an 
interpolation  function.  Inserting  eqs . (3)  and  (4)  into  eq. 

(2)  and  eliminating  stress  parameters  6^ by  using  stationary 
condition  for  the  functional  tt  with  respect  to  68^ , t can  be 
expressed  in  terms  of  generalized  displacements  q^  and  the 
stress  intensity  factor  . For  the  entire  structure,  terms 
associated  with  the  same  nodal  displacement  can  be  assembled. 
Again,  the  stationary  condition  of  n(q^,  k^)  with  respect  to 
6q^  and  6K^  leads  to  the  final  algebraic  equation  of  the  form, 
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Solving  the  above  equations,  the  generalized  displacements 
and  the  stress  intensity  factors  can  be  obtained.  It  is 
also  possible  to  consider  only  the  group  of  elements  at  the 
crack  front  and  form  a superelement.  The  eauations  are  of 
the  same  form  as  eq.  (5).  By  eliminating  {K^},  a stiffness 
matrix  [K]  of  the  superelement  given  by 


[K]  - [Krr]  - [Krs][KsJ1[Krs]T 


can  be  obtained.  Such  stiffness  matrix  can  be  used  in  the 
ordinary  matrix  displacement  method. 

Three  types  of  singular  element,  i.e.,  8-node,  12-node 
and  16-node  elements  are  constructed  (fig.  3).  The  16-node 
element  has  mid-side  nodes  on  all  edges  on  the  face  that 
intersects  the  crack  front,  whereas  the  12-node  element  has 
mid-side  nodes  only  on  edges  which  are  not  connected  with  the 
crack  front.  Table  1 lists  the  number  of  6's  used  in  the  var- 
ious elements.  Figure  3 also  indicates  typical  assemblages 
of  crack  elements.  The  12-node,  20-node  and  26-node  half 


elements  are  for  problems  which  are  symmetric  about  the  plane 
of  the  crack,  while  the  20-node,  36-node  and  46-node  super- 
elements are  for  general  asymmetric  problems. 

To  assess  the  accuracy  of  these  special  crack  front  ele- 
ments, a simple  test  has  been  performed.  Element  nodal  dis- 
placements are  calculated  from  the  known  analytical  solution 
that  corresponds  to  K ^ = K 2 = = 1.00.  Using  these  dis- 

placements as  input  data,  K-^,  K2  and  ate  then  computed 
from  the  numerically  generated  element  matrices.  The  re- 
sults obtained  by  using  different  types  of  singular  element 
and  by  using  different  half  elements  and  superelements  are 
listed  in  Table  1.  It  is  seen  that  they  all  compare  favorably 
with  the  correct  values.  Obviously  elements  with  more  nodes 
yield  more  accurate  solutions. 

The  elements  have  been  utilized  to  analyze  several 
fracture  test  specimens  of  typical  geometries,  i.e.,  a center 
crack  specimen,  a single  edge  crack  specimen  and  a double 
edge  crack  specimen  under  uniform  tension  along  the  faces 
parallel  to  the  crack  plane.  The  geometries  of  these  speci- 
mens and  mesh  subdivision  used  are  shown  in  fig.  (4)  and  (5) 
respectively.  By  the  double  or  triple  symmetry  of  geometry, 
only  one  quarter  or  one  eighth  of  these  specimens  are  con- 
sidered. Figures  (6),  (7)  and  (8)  show  the  variation  of  the 
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value  of  the  calculated  stress  intensity  factor  across 
the  plate  thickness.  Another  example  problem  that  has  been 
analyzed  with  this  special  element  is  a compact  tension 
specimen  shown  in  fig.  (9).  Because  of  the  double  symmetry 
of  geometry  only  one  quarter  of  the  specimen  is  needed  in  the 
analysis.  Figure  (10)  shows  the  mesh  pattern  used  for  this 
problem.  There  are  168  nodes,  85  elements  and  504  degrees  of 
freedom.  It  is  seen  that  the  pin  hole  is  removed  in  this  model 
and  the  loading  is  represented  by  a distributed  line  load 
normal  to  the  crack  surface. 

Figures  (11)  and  (12)  show  the  variation  of  the  stress 
intensity  factor  across  the  thickness  of  the  specimen 
for  two  different  geometries.  They  are  compared  with  the 
results  obtained  by  Yamamoto  and  Sumi  [6]  who  applied  a method 
involving  the  superposition  of  analytical  and  finite  element 
solutions.  The  problem  of  compact  tension  specimen  has 
also  been  analyzed  by  Barsoum  [7]  and  Tracey  [8]  using  isopara- 
metric block  elements  and  wedge  shape  elements  respectively. 
Both  elements  are  derived  by  assumed  displacements  which  con- 
tain the  singular  terms  for  and  k2  but  not  for  K^.  As 
indicated  in  Table  2,  by  using  a mesh  pattern  with  fewer 
degrees-of- freedom  than  that  by  Barsoum  and  Tracey,  the  pre- 
sent hybrid  stress  model  provides  results  which  are  closer  to 
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the  solutions  by  Yamamoto  and  Sumi.  The  latter  are  used 
here  as  reference  because  they  were  obtained  by  finer  meshes 
and  also  contain  a complete  treatment  of  the  singular  behavior 
at  the  crack  front. 

3.  Conclusions 

Preliminary  results  obtained  by  the  presently  developed 
special  3-D  crack  elements  are  satisfactory.  Further  studies 
are  needed  to  evaluate  the  merits  of  the  various  types  of 
elements  and  to  determine  the  applicable  ranges  for  the 
geometry  of  elements.  Finite  element  solutions  should  be 
obtained  for  problems  which  have  exact  solutions  to  be  used 
as  reference. 
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TABLE  1 

Evaluation  of  Different  Types  of  Element  and  Superelement 


Type  A 


No.  of  Nodes 


No.  of  8's 


K, 


8 


2 


.94250 


1.0325 


8 


24 


94241 


1.0287 


12 

16 

30 

45 

.98944 

1.0031 

1.0015 

.97760 

.99874 

.99988 

Type  B 


No.  of  Nodes 


No.  of  8's 


K, 


8 


24 


1.0126 


1.0279 


1.0019 


12 


5 


,99841 


99461 


99896 


1.0040 


1.0012 


1.0029 


16 


51 


1.0040 


1.0012 


1.0029 


Half  Element  (Type  A + Type  B) 


Total  No.  of  Nodes 

12 

20 

No.  of  Nodes 

8 (Type  1) 

8 (Type  2) 

12 (Type  1) 

12 (Type  2 ) 

No.  of  8's 

24 

24 

30 

45 

K1 

.98245 

.99608 

K2 

1.0136 

.99647 

K3 

1.0257 

.99886 

Table  1 (continued) 

Half  Element  (Type  A'  + B') 


Total  No.  of-  Nodes 

12 

No.  of  Nodes 

8 (Type  B 1 ) 

8 (Type  A' ) 

No.  of  B's 

24 

24 

.96286 

1.0497 

1.0257 


Complete  Super  Element  (Type  A + B + A'  + B') 

Total  No.  of  Nodes  20 

No.  of  Nodes  8 (Type  1,3)  8 (Type  2,4): 

No.  of  B's  24  24 

K,  .98025 


20 

36 

8 (Type  1,3) 

8 (Type  2,4) 

12 (Type  1,3) 

12 (Type  2,4) 

24 

24 

30 

45 

• 

98025 

.99312 

1. 

0290 

.99619 

J 

1.0257 


.99886 
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**K2[):  2-dimensional  plane  strain  stress  intensity  factor  solution 


FIG.  2 FINITE  ELEMENT  IDEALIZATION  AROUND 
CRACK  FRONT 


FAMILY  OF  THREE  DIMENSIONAL  HYBRID  STRESS  CRACK  ELEMENTS 


FRACTURE  TEST  SPECIMENS 


FINITE  ELEMENT  BREAKDOWN  FOR  THE  ANALYSIS  OF  FRACTURE 
TEST  SPECIMENS 


O 35  ELEMENT  MESH  SUBDIVISION  WITH 
20-NODE  SINGULAR  HALF  ELEMENT 

A 35  ELEMENT  MESH  SUBDIVISION  WITH 
12-NODE  SINGULAR  HALF  ELEMENT 

□ 85  ELEMENT  MESH  SUBDIVISION  WITH 

12-NODE  SINGULAR  HALF  ELEMENT 
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FIG.  7 SINGLE  EDGE  CRACK  SPECIMEN 
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FIG.  3 


DOUBLE  EDGE  CRACK  SPECIMEN 
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FIG.  9 


H = 0.6W 
h = 0.2  7 5W 
D = 0.2  5 W 
c = 0.2  5W 
2 t = 0.5  0W 


STANDARD  COMpACT  TENSION  SPECIMEN 
(ASTM  STANDARD  E-399-72) 


FINITE  ELEMENT  BREAKDOWN  FOR  THE  ANALYSIS  OF  STANDARD 
COMPACT  TENSION  SPECIMENS 


TWO-DIMENSIONAL  SOLUTION 


a YAMAMOTO ' S SOLUTION 

O 12-NODE  SINGULAR  HALF  ELEMENT 
COMPOSED  OF  TWO  8-NODE  BASIC 
SINGULAR  ELEMENTS 
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FIG.  .11  STRESS  INTENSITY  FACTOR  FOR  COMPACT 
TENSION  SPECIMEN 
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PIG.  12  STRESS  INTENSITY  FACTOR  FOR  COMPACT 
TENSION  SPECIMEN 


